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MTHSCOSH ; Floating Point Hyperbolic Cosine routi 16-SEP-1984 01:11:53 VAX/VMS Macro V04-00 Page 
1-006 gS Ei obs OF id i3a MTHRTL.SRCJMTHCOSH.MAR; 1 . 
09 1 -TITLE MTHSCOSH 3 Fhegeine Point Hyperbolic Cosine routine 
S «IDENT /1-006/ 3 Fite. MTHCOSH.MAR Edit: RNH1006 
8 5 ITITIIIITIIITITIIITI TLL LL TTILILILILTLi LITT Li TTT iLiiiiiiiriiiiiiTiittti tte 
;# * 
00 : 3* COPYRIGHT (c) 1978, 1980, 1982, 1984 BY * 
00 8 3* DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. * 
44 10 7 ALL RIGHTS RESERVED. * 
b+ ® 
00 11 ;* THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED * 
§ 1 3* ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE * 
0 13 ;* INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 
0000 14 ;* COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
0000 15 ;* OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
44 1 :* TRANSFERRED. x 
3 ® 
0000 18 ;* THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
0000 19 ;* AN OULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
44 y aN CORPORATION. * 
3 ® 
0000 § 3* DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS * 
0000 ;* SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. * 
0000 4 ;3* * 
0000 5 ;3* * 
0000 6 FERRARA RARER AREA AERA ARREARS 
0000 73: 
0000 Fe 
0000 9 ; FACILITY: MATH LIBRARY 
0000 Q ;++ 
a4 1 ; ABSTRACT: 
0000 : ; MTHSCOSH is a function which returns the floating point hyperbolic cosine 
0000 4 ; of its single precision floating point argument. The call is standard 
0000 2? : call-by-reference. 
0000 6 ; 
0000 7 ;- 
0000 8 ; 
0000 9 ; VERSION: 01 
8308 40 ; 
000 41 ; HISTORY: 
4! 5 3 AUTHOR: 
00 43; Peter Yuo, 29-Jun-77: Version 01 
Bp08 44 ; 
00 45 ; MODIFIED BY: 
44 46 : 
000 47 ; 
0000 48 ; 
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9 
80 


e 
urrent Edit History 6-SEP-19 MTHRTL.SRCJMTHCOSH.MAR; 1 
-SBTTL HISTORY ; Detailed Current Edit History 


; ALGORITHMIC DIFFERENCES FROM FP=11/C ROUTINE: none 
Edit History for Version 01 of MTHS$COSH 


0-2 MTHSSERROR Hie to MTHSS$SIGNAL. 
MTH$_... changed to MTH__.... 
Changed error handling mechanism. Put error result in RO before 
calling MTHSS$SIGNAL in order to allow user modify error result. 


0-3 Series changed from four term Taylor series to four term 

Chebyshev series, to improve accuracy. 18-May-1978; Bory Payne 

- — copyright notice and version number. JBS 16-NOV-7 

- Cha FLOOVEMAT to MTHSK_FLOOVEMAT. JBS 07-DEC-78 

- Remove SSRABDEF macro -- not needed. 495 Jenpace? 

- Add "*_"* to the PSECT directive. BS 21-DEC-7 

- Declare externals. SBL 17-May-1979 

- aeneee lower Limit for Chebyshev approximation from 2**-11 to 

- Eliminated second call to EXP for input values between .25 and 
12.5%ln2 by computing COSH(x) = (Z + 1/2)/2, with Z = EXP(ix!). 

- Et taloated second call to EXP for input values between 12.5*ln2 
an tln2. 

- Changed all final floating point divisions by 2 to interger 

subrtracts of 1 from the exponent field. 

Changed ok mask to excluded R5 - no longer needed. 

Extended maximum range from 87.69 to 128*\n2=88.72. 

zoaneges 2 °8r* for computing EXP(ixi-ln2) to reduce error. 


NOUS WN © OOO NA UE WWI OC OONOAUSE WO 
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-SBTTL DECLARATIONS ; Declarative Part of Module 


oe. 

wm 

aon 
— 
° 
we 
7 


INCLUDE FILES: 


: EXTERNAL SYMBOLS: 


: GBL 

sEXTRN MTHSEXP_RG 
-EXTRN MTHSK_FCOOVEMAT 
-EXTRN MTHSSSIGNAL 


0 3 

; EQUATED SYMBOLS: 
re g0k 03) 0 LF At OG. HI = MLE yn td 3; high 17 bit of ln2 + 2e*-17 
105CB 4 0 LF 606 LO = *x105CB7D0 : Ind - LF LOG gH 
O4F m4 0 LF-2 POWER_12.5 = “XO4F rr: : 2%212.5 = 579276188 
OF 3 4380 0 LF 1g7-L06.¢ = *X0F 334380 $ iF = Bb oe oets 
2174381 0 LF_128_L0G_2 = *X72174381 : 128%ln2 = 88.722839 
00004080 0 SF-1.0° = *F1.0 : 1.0 
00000004 value =4 3 value.rf.r 

0 

0 

0 


Soooo 


: PSECT DECLARATIONS: 


-PSECT _MTHSCODE PIC,SHR,LONG, EXE ,NOWRT 
3; program section for math routines 


So 
So 
S 
S 


OWN STORAGE: none 


CONSTANTS: 
COSHTAB: 
7381 3BB6 .WORD “0 35666,°0071601 : DECIMAL: 0.1391992E-02 
AABA £ A "WORD *0 7058 0125212 > DECIMAL: 0.416654E-01 
000 40600 WORD *0040006,0 > DECIMAL: 0. 3000000E +00 
000 4080 WORD  *0040200-0 : DECIMAL: 0.1000000E+01 


WWIWIAIR I RINPIDNINININININIDN 2 2 9 OP | “DOOD OOOOCOOOO00O0000000MMcmMccccDeD 


WN SO OONAOUES WN OOD NAULE WN 9 ODNAUE WN -2 0 OO NAME WN COONS WTP 
Sete . 
o- 
za 
oO 
n” 


aa a at ad td = = = = = 4 — — = 3 4) 2) — 2) od sd 2 8 2 ss 


SOOOCSSOSOOSOSOSOoOOOSoOOOoOOoOoOSoOoOSSoS 


00000004 COSHLEN = .- COSHTAB/4 


SOooooooooooooooooooooooo 
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i -SBTTL MTHSCOSH = Standard Single Precision Floating COSH 


FUNCTIONAL DESCRIPTION: 
COSH = single precision floating point function 
COSH(X) is computed as: 
< 2et-12 


3 If ix! , COSH(X) = 1, 
; IF Sait eC iki < 0342 Sag = Chebyshey series 
: If i ,' Bt S<,859 439 » set Z = EXP({X!), and compute 
0 : If gers <ixt < Bs. 029613 then COSH(X) = EXP({xX!)/2 
0 ; If 5-953 $i < iXt < 88.722839 then COSH(X) = EXP(iX!i=Ln2) 
8 ; If 2839 <= iX! overflow is signaled 
9 : CALLING SEQUENCE: 
0 ; COSH.wf.v = MTHSCOSH(x.rf.r) 
6 : INPUT PARAMETERS: 
00000004 0 LONG = 4 E getine Longuere multiplier 
00000004 . x = 1 * LONG : Contents of x is the argument 
0 IMPLICIT INPUTS: none 
6 ; OUTPUT PARAMETERS: 
VALUE: floating hyperbolic cosine of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 


SIDE EFFECTS: 


Signal: MTH$_FLOOVEMAT if 88.722839 <= {Xi with reserved operand in RO 

(copied to the signal mechanism vector CHFSL _MCH yet? by, LIBSSIGNAL). 

Associated message is: ‘FLOATING OVERFLOW IN-MATA LIBRARY’. Result is reserved 
operand -0.0 unless a user supplied (or any) error ee changes CHFS$L_MCH_RO/R1 


NOTE: This procedure disables floating point underflow, enables integer 
over fiow. 


ENTRY MTHSCOSH, “M<IV ,R2, R3, R4> 

: standard call-by-reference entry 

; disable DV (and FU), enable IV 
MTHSFLAG_JACKET :; flag that this is a jacket procedure in 


MOVAB G*MTHSSJACKET_HND, (FP) 
; set handler address to jacket 


ce cc ec ce ee ee ee a ee me ee ee ee ee a ee ce ec ee ee ee ee a ee eB ee ee ee a ed ed 
NOuw 


€00909 090909090909 SINISE A AA AAA UT EPP LEE 
DNA WN OOO NAME AN O OD NA UE WN 0 OO NAMES WO OONOULS WO 


DNIAMINININOOCOCOCOOOOCOCOCSOOOOOOOSOSCOOOOCCOOCOOCOOCOO OOOO OOOOOO OO OOOOO 


401C 
6D 00000000'GF 9 


nn —_ —_ 
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3; handler 
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case of an error in routine 
lf. teee convert signal to user PC 
and res 
50 4 BC 50 RO = X = @value(AP) 
0 2000 BF aA Sia 
F80 BF 0 861 compare {Xi with 0.25 
14 «18 branch if ix! >= 0.25 
3A80 BF 50 . 
0 
0 
0 


RO, #°XSF 
BGEQ  GEd_T0_0. 


Bicw2 od » 
CMPW 


; Xi < 0.25 
vant} BcEd GEG TO aMi> Branch if ix? se Seen] 
PIKE < Qeest2 
50 08 50 ; RO = 1.0 


MOVF S“#SF_1.0, RO ‘ 
T ; return with result = 1.0 


DD 
mm 
. 


: 2ee-12 =< IX! < 0.25 


GEQ_TO_2m12: 
MULF 


50 50 44 Q RO,RO ; Get ARG**2 for POLYF 
C4 AF s 7 OR he POLYF RO, #COSHLEN-1, COSHTAB ; RO = SUM(Ci®X*#*(2*{)) 
0 3; Last coefficient of 1.0 
8 3; provides eyerneng, 
04 8 RET 3; return with result in RO 
0 3 
3 0.25 =< ix! 
GEQ_T0_0.25: 
0F334380 8F 50 51 CMPF RO, #LF_127_L0G_2 3; compare {Xi with 127*ln2 
1c 14 BGTR GTR_THAN_127_L0G_2 : branch if iXi > 127*Ln2 


3 0.25 =< IK! =< 127#Ln2 


00000000" JSB FTHSEXP_R4 ; RO = EXP( 


BEF FAA WIAIAIG Ma RINPINIPINININININ 2 4 SS S| 2QOOOOOOCOCOCOCOWOOOOOOOOWOO 


WN SO ODNOAOUE WIN OOD NAUE WIN 9 OONA UE WI OO OO NOAU EWR OOONOUE W000 


F 16 004 ; Xd) 
04F346B85 8F § 51 004 CMPF =-RO, WLF-2 POWER_12.5 ; compare EXP(!X!) with 2*#12.5 
07 «14 «005 BGTR ONE_TERA_ONLY : If smaller EXP(={X!) is significant 
: : and must be added in 
51 08 50 4? DIVF3 RO, WSF 1.0, R1 > RI = EXP(=!x!) 
50 51 40 005 ADDF R1.R > RO = EXP(X) + EXP(-x) 
ne ONE_TERM_OWLY: 
50 0080 8F A2 005 SUBW © #*x0080, RO > RO = (EXP(X) + EXP(-x))/2 
04 $ RET : return with result in RO 
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44 
2a i 127ein2 < tx 
$ *in Ki 
ae 
6 49 GTR_THAN_127_L0G 2: 
72174381 8F 50 51 6 9 be CMPFO 6 #LF_128_L0G_2 ; Check for possible overflow 
19 «14 be 1 TR E : If greater, COSH(X) will overflow 
50 72804031 8F 42 26 SUBF #LF LOG_2 HI, RO ; RO = ixXi = (nigh order bits of Ln2) 
pp900 00'EF 19 2 MTHSEXP_ ; RO = EXPCIX!-LF_LOG_2_H1) 
51 50 1 8 DO 8F 4 7 4 MULF3 #LF_LOG_2_L0, RO, R1 : 
0 1 ae 0 2? SUBF R1, RO ; 
0 $ RET 
te 
4 58 
Bae 59 
084 60 ; 
084 61 ; 128*ln2 =< ixXi, error 
ak 
7E 00° 8F 9A 844 64 ERROR: MOVZBL #MTHSK_FLOOVEMAT, -(SP) ; condition value 
50 01 OF 78 0088 65 ASHL #15, #T, RO ; RO = result = reserved operand -0.0 
08C 96 3 goes to signal mechanism vector 
08C 6 ; CCHFSL_MCH_RO/R1) so error handler 
08C 68 ; can modify the resuit. 
00000000'GF 01 4=FB 008 69 CALLS #1, G*MTHSS$SIGNAL : gignal error and use real user's PC 
09 70 ; independent of CALL vs JS 
04 Si ¢] RET ; return - RO restored from CHF$L_MCH_RO/R1 
0094 7 
94 74 
0094 75 .END 


MTHSCOSH 
Symbol table 


FS 925 
COSHTAB 
OR 


Initialization 


{5g work ing 
bytes ( 
There were 1 


oe processing 
Syabol table sort 


Symbol table output 
Peect synopsis output 
Cross-reference output 
Assembler run totals 


et Limit was 900 
pages) of virtua 
pages of symbol table space allocated t 

source Lines were read in Pass 1, 

1 page of virtual memory was used to “define 1 macro. 


fererene 
fereeeee 
00000010 
eeerenee 
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= 00004080 
= 00000004 


Allocatio 


00000000 
00000094 
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5-3 0 NOPIC USR CON ABS LCL 
( 148.) O01 ¢ 1.) PIC USR CON REL 


$eoccccrcocececccoecoecon} 


LCL SHR EXE RD NOWRT NOVEC LONG 
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R 1 
R 1 
R 1 
a 1 
g 1 
x 1 
x 0 
% & 
x 00 
aa 01 
$¢ ee mreceeneuwaraceoce + 
: ;_Psect synopsis ! | 
PSECT No. Attributes | 
( 0°¢ 0.) NOSHR NOEXE NORD NOWRT NOVEC BYTE | 


BSR EEE 
1 Een Goes 
J GEES TERS 


producing 


memory were used to buffer the intermediate co 


hold 21 non-local and 8° local symbols. 
1 object records in Pass 2. 
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borer cen nena n anon cme somes} 


! Macro Library statistics ! 


tom nen SOseesso www ar sacoane a $ 


Macro Library name Macros defined 
_$255$DUA28:(SYSLIBJSTARLET.MLB;2 aT 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHCOSH/OBJ=OBJ$:MTHCOSH MSRC$:MTHJACKET/UPDATE=(ENHS :MTHJACKET) +MSRCS: 


AH-BT13A-SE EQUIPMENT CORPORATION 
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